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0^ ■ ABSTRACT 

P i 
■ ' We propose a formulation of an absorbing boundary for a quantum particle. The 

Ch ' formulation is based on a Feynman-type integral over trajectories that are confined by 

the absorbing boundary. Trajectories that reach the absorbing wall are instantaneously 

O^. terminated and their probability is discounted from the population of the surviving trajec- 

tories. This gives rise to a unidirectional absorption current at the boundary. We calculate 
the survival probability as a function of time. Several modes of absorption are derived 



Vh I from our formalism: total absorption, absorption that depends on energy levels, and ab- 

sorption of non-interacting particles. Several applications are given: the slit experiment 
with an absorbing screen and with absorbing lateral walls, and one dimensional particle 
between two absorbing walls. The survival probability of a particle between absorbing 
walls exhibits decay with beats. 



1. Introduction 

There are many different descriptions of absorption in electromagnetic tlieory, nuclear 
physics, solid state physics, and so on. These descriptions are mostly based on scattering 
theory. In scattering theory absorption occurs when a jet of particles hits an obstacle 
and is partially absorbed. The reduction in the population of particles is described in 
scattering theory by a complex potential. The introduction of a complex potential in 
scattering theory is supported by and large by phenomenological considerations |jl|. 

The introduction of a complex potential in scattering theory, that is, of phase shift, is 
analogous to the introduction of a complex dielectric coefficient in Maxwell's equations. 
There exist several models for the calculation of the complex part of the dielectric constant 
in Maxwell's theory. The basis for these models is the possibility to write equations 
of motion of a damped charged particle in an electromagnetic field and to identify the 
complex part of the dielectric constant as the damping coefficient [Q Q . Another approach 
to decay is based on the decoherent histories approach 0]. 

The purpose of this paper is to propose a Feynman-type integral to describe total 
or partial absorption of particles in a surface bounding a domain. The need for such a 
description stems from the reflecting property of potentials of all kinds, real and complex, 
flnite and inflnite. This property eliminates potentials as means of absorbing all particles 
that reach a given absorbing surface such as a photographic plate. Additional postulates 
are needed to incorporate absorbing surfaces into quantum mechanics. The method of 
complex potentials, as used in Maxwell's equations, does not carry over to quantum 
mechanics due to the fact that the wave function of a particle does not interact with 
the medium the way an electromagnetic wave does. In particular, in classical quantum 
theory, unlike in electromagnetic theory, the wave function does not transfer energy to 
the medium. 

In our approach to absorption, we postulate an absorbing boundary that absorbs 
instantaneously all Feynman trajectories at the instant they reach the boundary. The 
probability of these trajectories is discounted from the wave function of the surviving 
trajectories, much like the procedure adopted in the derivation of Fermi's Golden Rule. We 
flnd that the wave function of the surviving Feynman trajectories satisfles Schrodinger's 
equation with reflecting boundary conditions on the absorbing boundary. We deflne a 
discounted wave function whose total probability within the absorbing boundaries is the 
survival probability of the un-absorbed trajectories. We flnd the decay law of the survival 
probability for different modes of absorption: total absorption, absorption that depends 
on energy levels, and absorption of non-interacting particles. Total absorption is obtained 
if the wave function is discounted by the probability of the absorbed trajectories. In this 
mode of absorption the survival probability decays in time at a rate proportional to the 
total energy of the system with beats at Bohr frequencies. Absorption that depends on 
energy levels apportions different decay rates to modes corresponding to different energies 
in a such a away that the total discounted probability is the same as in the former case. 
Absorption of non-interacting particles produces a discounted wave function which is the 



product of the discounted wave functions of each particle separately. 

In Section 2, we formulate the postulates of absorption in a surface and explain the dis- 
counting procedure. In Section 3, we calculate the Feynman integral over the class of tra- 
jectories that are bounded within a given boundary and show that it satisfies Schrodinger's 
equation with total reflection on the boundary. In Section 4, we calculate the probability 
of Feynman trajectories that propagate into the absorbing boundary and obtain the sur- 
vival probability of the un-absorbed trajectories. In Section 5, we describe other modes of 
absorption. In Section 6, we calculate the absorption current at the absorbing boundary 
and consider the slit experiment with an absorbing screen. Section 7 contains the follow- 
ing examples: a particle between absorbing walls, a Gaussian wave packet incident on an 
absorbing wall, and the slit experiment with absorbing lateral boundaries. Experiments 
are proposed to examine the various decay modes. Finally, a discussion and summary are 
offered in Section 8. 

2. The postulates of absorption in a surface 

The two simplest types of instantaneous absorption at a wall are the absorption of 
all Feynman trajectories when they reach the wall for the first time, or absorption of 
trajectories that propagate across the wall. It was shown in [^] that Feynman integrals 
over the set of trajectories that terminate at a given wall produce wave functions that 
vanish on and beyond the wall. The resulting wave functions are continuous, though have 
discontinuous derivatives on the boundaries of their supports. It was shown in that for 
this type of initial wave functions the probability that propagates across the boundary of 
the support in a short time At is proportional to At^/^ so that a continuous discounting 
of the wave function by the probability of the trajectories that crossed into the absorbing 
domain beyond the wall leads to total reflection, much like in the Zeno effect |j^. It is 
shown below that the probability density that propagates into an absorbing wall in a short 
time At is proportional to At, so that discounting the probability of these trajectories 
leads to a decay law, consequently, we adopt the concept of instantaneous absorption at 
a wall for wave functions that vanish on and beyond the wall. 

More specifically, in our Feynman-type integral trajectories that propagate into the 
surface for the first time are considered to be absorbed instantaneously or reflected with 
a given probability which is a property of the interaction between the particle and the 
absorbing surface. The absorbed trajectories are therefore terminated at that surface. The 
population of the surviving trajectories is discounted by the probability of the absorbed 
trajectories at each time step. The instantaneous absorption rate is assumed proportional 
to the number (probability) of particles at the absorbing surface at a given moment of time. 
The proportionality constant is a characteristic length. A, determined by the absorbing 
material, e.g., photographic plate or fluorescent screen, and by the absorbed particle, e.g., 
electron or neutron. The length A is assumed to depend weakly on the absorbed particle's 
energy so, in this paper, we approximate it by a constant. Under this assumption A can be 
determined experimentally by measuring the absorption rate of a single energy level and 



than be used to calculate the absorption rate at any other energy level, in some energy 
range, or the simultaneous absorption of several levels in this range. The constant A can 
be assumed to contain a multiplicative factor that represents the probability of absorption 
of a trajectory when it hits the wall. 

In the context of the decay of an energy level into a continuum, usually described by 
the Wigner-Weiskopf theory and Fermi's Golden Rule ||^, the characteristic length A in 
our theory, can be related to the size of the support of the perturbing Hamiltonian. In 
the context of a-decay in nuclear physics, the length A can be related to the width of the 
potential barrier p. 

A detailed theory of A can be expected to involve the atomic detail of the absorbing wall 
and the details of the interaction between the particle and the atoms of the absorber. In 
contrast, our theory involves the trajectories of the particles before they hit the absorbing 
wall so that A is an external input into our theory. This situation is analogous to the 
theory of frequency dependent refraction coefficient in electromagnetic theory where the 
frequency dependence is determined from a detailed analysis of the interaction between 
the electromagnetic wave and the atoms of the medium. Thus the frequency dependent 
refraction coefficient is an external input into the Helmholtz equation. In this theory, 
the knowledge of the response to a single frequency does not determine in a simple way 
the response to a wave containing several frequencies. Only in the case of a frequency 
independent dielectric constant its measurement at one frequency allows to predict the 
response to a wave containing several frequencies. In different theories of absorption the 
parameter A can depend on any number of physical properties of the system, such as the 
energy of the system or some average energy per mode, and so on. 

In general, if the trajectories are partitioned into two subsets, the part of the wave 
function obtained from the Feynman integral over one subset cannot be used to calculate 
the probability of this subset, due to interference between the wave functions of the two 
subsets. However, in the physical situation under consideration, such a calculation may 
be justified as follows. Our procedure, in effect, assumes a partition of all the possible tra- 
jectories at any given time interval [t, t+ At] into two classes. One is a class of restricted 
trajectories that have not reached the absorbing surface by time t + At and remain in the 
domain, and the other is a class of trajectories that arrived at the surface for the first time 
in the interval [t,t+ At]. We assume that the part of the wave function obtained from 
the Feynman integral over trajectories that reached the surface in this time interval no 
longer interferes with the part of the wave function obtained from the Feynman integral 
over the class of restricted trajectories in a significant way. That is, the interference is 
terminated at this point so that the general population of trajectories can be discounted 
by the probability of the terminated trajectories. This assumption makes it possible to 
calculate separately the probability of the absorbed trajectories in the time interval [t, t+ 
At]. 

The assumptions discussed above can be summarized in the form of the following pos- 
tulates: 



I. The absorption of Feynman trajectories represents absorption of actual trajectories of 
particles. 

II. The population of Feynman trajectories can be discounted by the probability of the ab- 
sorbed trajectories. 

3. The Feynman integral over bounded trajectories 

If a quantum particle is constrained in space to a finite (or semi finite) domain, 
the Feynman integral has to be confined to Feynman trajectories that stay forever in this 
domain. This implies the following modification in the definition of the Feynman integral. 
The function space is now the class 

a^b = {x{-) eC[0,t] I a < x (r) < 6, < r < t} 

and the definition of the Feynman integral over the class aa^b is 

Ka,b{x,t) = ^ expl^^S[x{-),t]jVx{-) 

= lim UN ■■■ exp <^ tS{xq, ...,XN,t)\ J| dxj. (3.1) 



Next, following the method of [T^, we show that Kafi{x, t) satisfies the Schrodinger equa- 
tion and determine the boundary conditions at the endpoints of the interval [a, 6]. We 
begin with a derivation of a recursion relation that defines K{x,t). We set 
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fb fb r i ] "'"^ 

KN{xN,t) = aNj ■■■] exp|-5(a;o,...,XAr,t)| n c?Xj, (3.2) 

then by the definition ( p.l|) , K{x,t) = \im.j^_,ooKj^{x,t). The definition (|3.2|) implies the 
recursion relation 

KNix,t)= (3.3) 

J I / ---J ^ ±2^-V(.^Af \KN-l{xN-l^N-l)dXN-l. 



I 2mh/\t J J a ^^^ 1 h 



2At ^ ' 



The following derivation is formal. A strict derivation can be constructed along the lines 
of [|10|. We expand the function Kjq_i{xN-i,tN-i) in (|3.3|) in Taylor's series about x to 



obtain 



KM{x,t)= (3.4) 
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K (r f ] (r T -. ^KN-iix^tN-i) 
J\N_l[X, iN-l) — [X — Xn-i) ^ h 



\2 



1 29^N-lix,tN^l) 

2(^-^^-) a^5 



+ {{X - XN-lf) 



dx 
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The integrals in eg. ( ^.41 ) are evaluated separately for x inside the interval [a,b] and on 
its boundaries. When x is inside the interval, the integrals become the Fresnel integrals 
over the entire line in the limit At ^ 0. This recovers the Schrodinger equation inside 
the interval. When x = a,b, the first integral in eq. ( |3.4| ) becomes in the limit At ^ the 
Fresnel integral over half the real line while the other integrals vanish. We obtain 

K{a,t) = ^K{a,t), (3.5) 

hence K{a,t) = and similarly, K{b,t) = 0. More specifically, consider the normalized 
n-th moment of the Gaussian integral (n = 0, 1, 2) 

a"m„(a, a, 6, y) = < ^ > / {x — y)" expl—ia{x — yy\ dx, 

where 

m 

2;iAt 

We change variable to 

y/a{x — y) = u (3.6) 

and get 

r 1 ] 1/2 fV^ib-y) f -, 

mn{<y,a,b,y) = < — > / u'^exp<—iu > du. 

Note that the limits of integration become 

\fQ.{a — y) — > — oo, \fa{b — y) ^ oo as a ^ oo 
so that in the limit At ^ 0, that is, as a — > cxo and for a < y < b, we obtain 

r^/a{b~y) , 

mn(a,b,y) = lim m„(a, a, 6, w) = / u^exp\—iu > du 



r u^'expl-iu^} du. (3.7) 

J — CO ^ ^ 



Note that mi = while mo = m,2 = 1- Now, transferring K^-i to the left hand side 
of equation ( |3.4|) , dividing by At, taking the limit, and applying eg. ([3171) , we recover the 
Schrodinger equation. 

At the boundary point y = b, the substitution ( p.6|) transforms the upper limit of 
integration to 0. In the limit At -^ 0, we obtain 



m.n(a, b, b) = lim r7i„(a, a, b, b) = lim / -u" exp I —iu^ \ du 

= / u"' exp l—iu^\ du. (3. 



That is, the Fresnel integrals are extended only over half the line. Thus, we obtain that 
the coefficient of Kat-i in eq.( p.4|) is 

mo(a,6,6) = -, 

while that of the first derivative of Kn_i is 

a~^m.i{a, a, b,b) ^ 
and that of the second derivative of K^^i is 

a~^r7i2(«, CL, b, b) — > 0. 



Thus, setting a; = 6 in eq. (|3.4|) and taking the limit At -^ 0, we obtain eq.( p.5|) , which in 
turn implies the boundary condition eq.( |3.10| ). 



When a < X < b, the derivation given in [|l^ leads to the Schrodinger equation. Thus 



th ^-^-^ = ^\ ' +V{x)K{x,t) for a < X < 6 (3.9) 

ot 2m, ox^ 

K{a,t) = K{b,t) = fort>0 (3.10) 

K{x,0) = 6{x-xi) foia<x<b. (3.11) 

In the case a = — oo, b = oo every x is an internal point so that the Schrodinger 

equation is satisfied for all x and thus the Feynman integral ( p.l| ) is equivalent to the 
Schrodinger equation (|3.9| ) and the initial condition ( p.ll| ) on the entire real line. 



This result shows that bounded trajectories imply in effect an infinite potential barrier 
on the boundary. In contrast, the Wiener integral over bounded trajectories leads to the 
diffusion equation with zero boundary condition and to the decay of the population. 
This means that the boundary is absorbing. This difference between the Wiener and the 
Feynman integrals over the same class of functions is illustrative of the different roles that 
trajectories play in quantum and classical theories. 



The same result has been derived in |]TT| and |[T^. In [|ri[] trajectories are extended 
into the entire hne as periodic functions. The analysis of this section provides a simpler 



derivation of the result. The derivation in [12| is done by converting the Wiener integral 



into the Feynman integral by analytic continuation of t into the imaginary axis. 
4. Feynman integrals with absorbing boundaries 

Assume now that a trajectory that reaches the boundary x = a ot x = b for the first 
time is instantaneously absorbed. This means that the wave function outside the interval 
(a, b) vanishes identically and the population inside the interval (a, b) is reduced at a 
rate determined by the current at the boundary, as described below. The vanishing wave 
function outside the interval (a, b) expresses the assumption that once outside the interval 
the particle no longer participates in the quantum evolution of the particles inside the 
interval. This may occur, for example, in the scattering of particles on a target [e.g., a 
nucleus). Particles absorbed in the nucleus are discounted from the scattered population. 
Another example is that of a particle that enters a bath, such as a photographic plate, 
and leaves an irreversible trace. Consequently, its quantum interaction with the particles 
inside the interval becomes negligible. Also in this case the population inside the interval 
is discounted at the rate particles are absorbed. 

The absorption process at a given time t is described as the limiting process as At ^ 
of the propagation of a trajectory that survived in the interval [a, b] till time t to a 
boundary point in the time interval [t, t + At]. In order to incorporate this behavior into 
the Feynman formulation, we adopt the procedure that leads to the Feynman-Kac formula 
for the probability density function for a diffusion process with a killing (absorption) 
measure (see, e.g., [0). We consider separately the trajectories that reach the boundary 
in the time interval [0, At] , then those that survived till At, but reach the boundary in the 
time interval [At, 2At], and so on. In each time step, the total population of trajectories 
has to be discounted by the probability of the absorbed trajectories. This leads to a 
modified expression for the discretized Feynman integral. 

First, we calculate the discretized Feynman integral in the time interval [0, At], 

?/'i(x. At) = 1^— ^— I J iJo{xo)expl-S{xo,x,At)\ dxQ. 

Therefore, the probability density of finding a trajectory at x = a in the time interval 
[0, At] is |?/'(a, At)| , and there is an analogous expression for the probability density of 
finding a trajectory at x = 6 in the time interval [0, At]. 

For simplicity, we assume that b = and consider the interval [—a, 0], where a > 0. In 
any time interval [t, t + At], the probability density propagated from this interval into the 
absorbing boundary at x = is calculated next. We begin with an initial wave function 
ip{x,t) that is a polynomial 



N 



g(x,t) = ^g,(t) 



X 



J 



in the interval [— a,0], such that Q{—a,t) = Q{0,t) = and ilj[x,t) = otherwise. The 
free propagation from the interval [—a, 0] is given by 



iij{y,t + At) = ^-^ 



Q[x, t) exp < — _. . — } ax. 



2'7iihAt J~a [ 2hAt 

The boundary condition at the left end of the support of ilj{x,t) is written explicitly as 



N 






Setting 



a 



hAt 



m 



the probability mass propagated freely into the absorbing boundary point x = in time 
At is given by 



l*'"-' + ^"'-2™ 


J-a 


2 


Dy setting x = ^/a^ to get 


|^(0,t + At)|^ = i- 


/: 





(4.1) 



First, we evaluate the inner integral, 



^ /-O 



All limits of the type 



lim «-(^+i)/2 / x^e'-V^-dx, 



are understood in the sense 



lim «~(^+^)/' r x^e'"'/2"dx= lim lim a"(^+^)/2 Z" x^e^''^'^'''^^'' dx. (4.2) 

'^0+ J~a^ e->0+Q->0+ J-a^ 



The first term in the sum Jtv gives 

rO 



ia'^/2a 



-a/y/a ^ 

The second term gives 



Setting 



J~a/Ja 



integration by parts gives the recursion relation 

Proceeding by induction, we find that for j > 2 

Sj = O (v^e*"'/2a 

Now, using the definition ([4 .21), we find that 



2imi|^(0,t + At)|^ = i-|gi(t)r. 



(4.3) 



Keeping in mind that 



giW 



OHM 
dx 



we can write 



|^(0,t + At)| 



2 ^At 



2Trm 



dij{0,t) 



dx 



fnAt\ ^ 

+ o for a < 1. 

\2nmJ 



(4.4) 



Returning to the propagation from an interval [a, b] to its absorbing boundaries, the 
Feynman trajectories in the time interval [0, At] consist of those that have not reached 
the absorbing boundaries and of those that have. We calculate the Feynman integral 
separately on each one of these two classes of trajectories. The integral over the first class 
is the one calculated in Section 3 above. That over the latter class is the integral calculated 
in the above paragraphs in this section. According to our assumptions, trajectories that 
propagate into the absorbing boundary never return into the interval [a, b] so that the 
Feynman integral over these trajectories is supported outside the interval. On the other 
hand, the Feynman integral over the bounded trajectories in the interval is supported 
inside the interval. Thus the two integrals are orthogonal and give rise to no interference. 

The Feynman integral over the bounded trajectories represents the wave function 
conditioned on not exiting the interval in [0, At]. According to eq. (|4.4| ) and to our as- 
sumptions, the probability defined by the Feynman integral over the trajectories that 
propagated into the absorbing boundary in this time interval is given by 



Pi(At) = A(a)|V^i(a,At)|' + A(6)|^i(6,At)|Vo(At) 



^At 
27rm 



A(a) 



g^i(a,0) 
dx 



+ m 



dx 



o(At) 



where A (a) and A (6) are the characteristic lengths at x = a and x = b, respectively. The 
survival probability in this time interval is 



Si{At) = l-Pi(At). 



(4.5) 



Next, we calculate the discretized Feynman integral in the time interval [At, 2At] and 
again break it into the same classes as above. The discretized integral over the class of 
bounded trajectories in the interval [a, b] is 



ip2 (x, 2 At) 



m 



.] 



1/2 



2'KihAt} 



/ '?/'i(xi. At) exp <^ —S{xi, X, At) > dxi. 



The probability calculated from the Feynman integral over the absorbed trajectory in the 
time interval [At, 2At] is 



P2(2At) = A(a) |^2(a, 2At)|' + \{b) |^2(6, 2At)|' 



o(At). 



This is the conditional probability of trajectories that propagate into the absorbing bound- 
aries in the time interval [At,2At], given that they did not reach the boundaries in the 
previous time interval. Thus the discretized survival probability in the time interval 

[0, 2 At] is 



82(2 At) 



Pi(At))(l-P2(2At)) 



Proceeding this way, we find that the discretized Feynman integral over the class of 
bounded trajectories in the time interval [0, A^At] is 



^pNix,NAt) 

m 1 1/2 rh 

27Ti 



(4.6) 



77, ^ 1/2 /-fc r 2 ^ 

— -| J ijN^i{xN-i,{N -l)At)expl-S{xN-i,x,At)\ dxN-i- 



As in eq. (|3.3|) , we find that iPn{x, NAt) -^ ip{x, t) a.s N ^ 00, where ip{x, t) is the solution 
of Schrodinger's equation in (a, 6) with the boundary conditions ilj{a,t) = ip{b,t) = 0. 

The probability that propagates into the absorbing walls in the time interval [{j — 
l)At, jAt] is given by 



P,-(jAt) = A(a) \^,{a,jAt)f + X{b) \^j{b,jAt)f + o{At). 
It follows that the survival probability of trajectories inside the interval is 



N 



^(t)= hm Y[il-P,ijAt)). 



N^oo 



(4.7) 



(4i 



i=i 



According to eqs.([4.4|) and ( |4.71 ), Pj{jAt) is given by 



P,(jAt) 



hAt 

211171 



A(a) 



d 

— ^,_i (a, (j - l)t) 



+ m 



d_ 
dx 



7A,_1 (6, (j - l)t) 



0(1) 



10 



so that eq.( [4.8|) gives the survival probabihty 



S{t) = exp 



h 



TCTTl Jo 



X{a) 



§-J{a^') 



m 



i.^^'^'' 



dt' 



(4.9) 



The wave function of the trajectories that have not been absorbed by time t is the 
wave function of a particle, conditioned on not reaching the absorbing boundary by time 
t. The conditioning renormalizes the wave function inside the domain at all times and 
thus it remains tp {x, t). The survival probability is 1 — P{t). If a particle is known not to 
have been absorbed by time ti, its survival probability till time ^2 > ^i, denoted S{t2,ti), 
is given by 



S{t2,ti) = exp 



h 



Tvm Jt 



t2 



\{a) 






+ m 



|^(M) 



dt 



(4.10) 



The survival probability S{t) and the wave function il){x,t) can be combined into a 
discounted wave function 



■^{x,t) = ^S{t)ilj{x,t). 



(4.11) 



The discounted wave function should be used for all purposes if an absorbing surface is 
present in the system. It is normalized to S{t) and decays in time. The the wave function 
ip{x,t) is recovered from \E'(a;, t) by conditioning on survival of the particle by time t. 

5. Other modes of absorption 

There are different phenomenological descriptions of absorption and decay in the lit- 
erature. These may include different absorption rates at different energies, absorption 
of two non-interacting particles at a wall, and so on. It is not a priori obvious which 
mode of absorption applies in a given physical situation without a detailed analysis of the 
absorption mechanism. We show below that some of the phenomenological descriptions 
of absorption can be derived from our formalism. 

5.1 Absorption in energy windows 

Some of the descriptions are based on the premise that each mode is absorbed (or 
decays) at a rate proportional to its energy. This type of absorption can be obtained 
from the model of absorption of Feynman trajectories at an absorbing wall. Since the 
trajectories that reach the wall in time At cannot be separated into those that belong 
to one mode or another, the overall discounted probability has to be a function of the 
rate and probability of propagation of trajectories into the wall and the overall energy of 
the system. The way the discounting is spread among the different components of the 
surviving wave function can be chosen according to any number of criteria. 



11 



Therefore, we adopt the following absorption principles: the overall probability dis- 
counted at each time step is a function of the rate of propagation of probability into the 
wall and the average energy of the system; the probability discounted from each mode 
is, in addition, a function of the energy of that mode. According to these principles, the 
factor A is allowed to be a function of the propagation rate, that is. 



A = A 



h 



Tim 



d_ 
dx 



^(0,t) 



In addition, it can depend on the average energy, {£), defined as follows. Consider a single 
absorbing wall placed at the origin and assume that the wave function in the presence of 
an infinite potential wall at the origin is 



ij (x, t) = Y^ ttk^pk (x) exp I --rEkt > 



where the eigenfunctions ipk (x), corresponding to energies Ek, form a complete set. The 
sum is replaced with an integral if the spectrum is continuous. The normalization condi- 
tion is 

oo 

El |2 1 
\ak\ = 1. 



fc=i 
The energy of a mode is defined as 

^fe = l^fcl Ek 
and we define the average energy per mode as 



(5.1) 



k=l k=l 



Thus the factor A has the form 

A = A 



h 

Tcm 



il**"-" 



,{£)]■ 



The overall probability propagating into the wall in the time interval [t, t + At] is given 
by 



P{t) 



Xh 



TTTTl 



d_ 
dx 



V^(o,t) 



At + O f At^/^l . 



(5.2) 



Now, we spread this discounted probability among the different modes. The mode k 
is discounted by 



flfcV'fc {x) exp \ -^Ekt M 1 - IT 

' h } \ 2 7im 



d_ 
dx 



^(0,t) 



At + O f At^/^ 
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where A^ is defined below. The value of Xk has to be chosen so that the overall discounted 
probability is as given in eq. (|5.2| ). The total surviving population is 



/ Yl ^k-ipk (x) exp I -T^kt I 1 



2 nm 



§-jm 



At + (At^^/^ 



dx 



2^ \ak\ 1 - Afc 



fc=i 



Tcm 



l^*M 



At + (At^/^ 



so that the overall discounted probability is 

2 



El |2 \ "■ 



k=l 



Tim 



li>(o,t) 



\h 



At + (At'/') = — 



TCTTl 



^:-"(0A 



At + (At^/2^ . 



Simplification leads to 

We set 

so that eq.( |573D becomes 



/ , \o-k\ Afc — A. 



(5.3) 



k=l 



Xi. = Xai 



El |2 1 

|afc| ak = 1. 

A:=l 



To obtain absorption rates proportional to energies, we have to choose 

^k 



ttfc 



It follows that 



Afc — A 



^k 



{^) 



(5.4) 



and the absorption law for each mode reduces to 



exp 
Now, we choose A so that 



A 



{S) Tim 
d 



!*■ '"•«'' 



dt'\. 



h 



TTm 



dx 



iPiO^t) 






we obtain the discounted wave function of the surviving trajectories in the form 
^ (x, t) = Y^ akipk (x) exp I -T^kt I exp I --^t i . 



(5.5) 
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The survival probability of mode k is 



Sk{t) = |afc|^ exp|-ytl 
and the total survival probability is the sum 

Ssit) = I" |^(x,t)|' dx = f; \ak\^ exp |-^t| . (5.6) 



fc=i 



This mode of decay is similar to that in diffusion, where each mode decays separately at 
a different decay rate. 

Now, we consider an interface between two media, one of which absorbs particles with 
energies in a given band El < E < Ejj and reflects all particles with energies outside 
the band. This situation is described by an infinite potential at the interface for particles 
with energies outside the band and an absorbing wall for particles with energies inside 
the band. 

According to the formalism developed above, all particles see effectively an infinite 
wall at the interface, but only those with energies in the band are absorbed. To describe 
the absorption process, we separate the initial wave function iPq{x) into its projection, 
ipband{.x) , ou the subspacc spanned by eigenfunctions corresponding to the infinite wall 
whose eigenvalues (energies) are in the band, and to its projection, '?/'c(a;), on the orthog- 
onal complement of this subspace. The wave function, ip{x,t), is also decomposed into 
analogous projections, 

Ipix, t) = i^bandix, t) + ?/'c(x, t), 

where ipband{x,t) decays according to the above absorption formahsm, whereas ipc{x,t) 
evolves according to the Schrodinger equation with zero boundary conditions on the in- 
terface. The Afc for this case are chosen by eq.( [5.4| ) for energies in the band and outside 
the band. The mean energy per mode is defined only for the modes that are absorbed, 
that is, 

fceband 
5.2 Absorption of two non-interacting particles at a wall 

The Hamiltonian for two non-interacting particles has the structure 

H{xi,Pi;X2,P2) = Hi{xi,pi) + H2ix2,P2) 

so that the wave function has the form 

^{Xi,X2,t) =ijJi{xi,t)^2{x2,'t)- 

If an absorbing wall for both particles is placed at x = and initially both particles are 
to the left of the wall, the configuration space is the third quadrant in the (xi, a;2)-plane. 
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The absorbing boundary, denoted dD, consists of the negative xi and X2 axes. Each 
particle can have a different A so that a two dimensional analog of the calculation of the 
absorption rate gives the instantaneous rate 



7rmi 



5V^i(o,t) 



dxi 



\ll)2{x2,t)\ dX2 



\,n /-o 



vrm2 



5V^2(0,t) 



dxi 



\i)i{xi,t)\^ dxj 



(5.7) 



TTfrii 



dMO,t) 



dxi 



+ 



\2Tl 

vrm2 



9^2 (0,t) 



dxo 



Thus the discounted wave function is given by 

<^{Xi,X2,t) = ^i{Xi,t)^2{x2,t), 

where 



\E'j {xj , t) = ipj {xj , t) exp 



Ao/i 



7rm,- 



9^,(0,t) 



dxj 



(J = l,2) 



Thus, the population decays together at a rate that is the sum of the two rates. The 
decay law is different than that of a single particle in a superposition of two states in that 
there are no beats in the decay law of two independent particles. 

It is apparent from the above examples that different modes of decay can be obtained 
from the model of an absorbing wall. It is not a-priory clear which mode of decay is 
appropriate for a given physical system. The choice of the decay mode depends on the 
particular physical system. 

6. The absorption current and the sht experiment 

An absorbing boundary engenders an absorption probability current into the boundary 
0. The simplest definition of the absorption current at a single absorbing point, x = 0, 
in one dimension, is (see eq. (^^) ) 



J{0,t) 



d 



dt 



[1 - S{t)] 



(6.1) 



a^(o, t) 



dx 



exp 



\h /■* 



vrrri jo 



9V^(0,t') 



dx 



dt' 



In higher dimensions the discounted wave function in a domain V in the presence of 
an absorbing boundary F is given by 



^(x, t) = ?/^(x,t)exp 



\n 



27rm Jo Jt 



(9^(x', t') 



dS^i dt' 



(6.2) 
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where ip{^, t) is the solution of Schrodinger's equation in T> with zero boundary condition 
on r and n is the unit outer normal to the boundary. Adopting the interpretation of the 
squared modulus of the wave function as the probability density of finding a particle at 
the point x at time t (whatever that means), the discounted wave function can be used 
to calculate the joint probability density of surviving by time t and finding the particle at 
X at the same time. The squared modulus of the wave function, conditioned on surviving 
by time t is found by dividing the joint probability density |\l/(x, t)| by the probability 
of the condition, S{t). In the multi-dimensional case at hand 



S{t) = exp 



Xh 



■Km Jo 



a^(x',t') 



dn 



dS-x' dt' 



(6.3) 



It follows from eqs.( |6.2| ) and ( |6.3| ) that the conditioned wave function is ip{:x.,t). 



Generalizing the definition (|6.1| ) to higher dimensions, the total current at the absorb- 
ing boundary is defined as 



Mt) = |[i-5W] 



Xh 

7im Jr 



9^(x,t) 



dn 



dS:ii' exp 



Xh /•< 



TCTTi Jo Jr 



dn 



(6.4) 



dip{x',t') 



dS-x' dt' } . 



The normal component of the multi-dimensional probability current density at any 
point X on r is given by 



^(x,t)-n(x)|r = 

Tim 



dH^,t) 



dn 



Xh /•* 



exp 



TCTTi Jo Jr 



dip{x',t') 



dn 



rfSx' dt' } . (6.5) 



Next, we consider the slit experiment with an absorbing screen {e.g., photographic 
plate or a fluorescent screen). In this case, unlike the usual assumption in quantum 
mechanics W^, the pattern that appears on the absorbing screen cannot be the squared 
modulus of the wave function in the entire space (with or without a finite potential), as the 
above analysis implies. It cannot be the modulus of the wave function with an absorbing 
boundary because the wave function vanishes on an absorbing boundary. According to 
the absorption principles, the pattern on the screen represents the probability density 
that propagates into the screen. This probability density is the unidirectional absorption 



current on the screen. Thus, the instantaneous pattern at time t is J7'(x, t) 
the cumulative pattern, after the arrival of many particles, is 



n x 



and 



Ji^) 



n X 



J{^,t) 



n X 



dt. 



(6.6) 



The instantaneous pattern i7'(x, t) ■ n (x)|p is obtained as the histogram of points on the 
screen collected at time t after releasing the particle in the slit. The cumulative current 
density J7'(x) ■ n (x)|p is the histogram of points on the screen collected at all times. 
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We consider the following experimental setup for the slit experiment with an absorbing 

screen. A planar screen is placed in the plane a; = and another screen it placed in the 

plane x = xo and it is slit along a line parallel to the 2;— axis. Due to the invariance of 

the geometry of the problem in z the mathematical description of the slit is, following 

1^ , an initial truncated Gaussian wave packet in the (x,?/)— plane, concentrated around 



the initial point, (xo,0). To describe the interference pattern on the screen, we assume it 
is an absorbing line on the y-axis and apply the formalism developed above. The wave 



function, as given by our formalism, evolves from the initial packet according to eq.(3.2 
as 



'^ix,y,t) =ip{x,y,t)exp 



Xh 



27im 



di'{0,y',t') 



dx 



dy' dt' 



where ilj{x,y,t) is the solution of the Schrodinger equation in the half plane a; > with 
■0(0, y,t) = and ip{x, y, 0) is the given initial packet. 

The probability current density on the screen is given by 



J{0,y,t) = 

rem 



dij{0,y,t) 



dx 



t roD 



exp 



Xh 



Tim Jo J~oo 



dij{0,y\t') 



dx 



dy' dt' 



(6.7) 



This is the probability density of a collapse of the wave function occurring at the point y 
on the screen at time t. The total current 



Jiy) = — / 

TTm Jo 



dtp{0,y,t) 



dx 



exp 



Xh 



vrm JO J-00 



dtp{0,y',t') 



dx 



dy' dt' > dt (6.t 



is the probability density that the collapse of the wave function occurs at the point y on 
the screen (ever). If the initial distribution of velocities in the x direction is concentrated 
about vq, the density J^{y) is approximately the same as J'{0,y,t), where Jll 



In a real experiment the measurement is neither instantaneous nor infinite in time. 
That is, an integral over a finite time interval is observed rather than (|6.7| ) or ( |6.8| ). If a 
packet of particles is sent out eq.( |6.7|) is the probability density of Feynman trajectories 
that propagate instantaneously at time t into the point (0, y) in the screen and is seen as 
the density of light intensity on the (ideal) fluorescent screen at time t. The function (|6.8|) 
represents the cumulative (in time) probability current density of Feynman trajectories 
absorbed in the wall and is seen as the density of the trace the initial packet eventually 
leaves on the screen {e.g., photographic plate). 

To determine the patterns (|6.7| ) and ( |6.8| ), we have to calculate first the two-dimensional 
wave function with zero boundary condition on the y-axis. It can be written as 

0(a;,2/,t) =0^(x,t)V'^(|/,t), 
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where ip^{x,t) is the function K{x,t) in eqs.(p.9|, |3.10| , p.ll| ) gaven above and by eq.( |7.4| ) 
geven below. 

Note that according to our theory ip^i^O, t) = 0, because the absorbing screen is placed 
at X = 0. The component ip'^{y, t) of the wave function is given by the freely propagating 



Gaussian slit |T4| 



r{y,t) 



Evaluation of the integral gives 



oo g-z^/2cr2 ^-im{y-z)'^/2ht 



2'Kia. 



y J2niht/m 



dz. 



^\y.t) 



2+2 



1 / rt 

27rcr„ I a'},m? ^ , 



-1/2 



exp < 



2+2 



UH 



aim- 



+ < 



If (Ja; << |xo|, the upper limit of integration in eg. (^^) can be replaced by cxd with a 
transcendentally small error. 

According to eq. (|6l7| ), the instantaneous absorption rate at time t at a point (0,?/) on 
the screen is given by 



J{^,y,t) 



dK{0,t) 



dx 



r{y,t) 



(6.9) 



where \dK{0, t)/dx\ is given in eg. (|7. 5|) below. Thus the pattern on the screen at time t 



is given by J'{0,y,t) in eg. (|6.9| ). 

To compare eg. (|6.9|) with that given in ||I^, we reproduce the derivation of [|1^ with 



an initial two-dimensional Gaussian wave packet. The result gives the wave function as 

iPf{x, y, t) = ipp{x, t)ilj^{y, t) 
and probability density at the screen at time t as 



\Mo,y,t)f 



277(71, \ aim? ^ , 



-1/2 



exp 



Xn 






i^'iy,t) 



(6.10) 



Comparing eguation ( |6.10|) with eg.( |6.9|) , we see that the pattern predicted by Feynman's 
theory differs from that in the present theory by a time dependent factor only. The 
instantaneous intensity of the diffraction pattern in the absence of an absorbing screen, 



given in [III], is defined as \ipF{0,y,t)\ . Thus, the introduction of an absorbing screen. 



according to these interpretations, gives the relative brightness as 



A 



dip{0,y,t) 



dx 



\i^F{0,y,t)\ 



2 ' 
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which is a function of time only. The decay in time of the quotient reflects the fact that 
the absorbing screen depresses the entire wave function in time. Thus, the part of the 
packet that arrives later is already attenuated by the preceding absorption, relative to the 
unattenuated wave function in the absence of absorption. 

7. Examples 

In the following examples possible experiments are discussed, that can be used to ver- 
ify the predictions of the above theory in the different modes of absorption. The examples 
include a particle between absorbing walls, a free particle incident on an absorbing wall, 
and the slit experiment bounded by absorbing walls. The latter indicates that the pattern 
on the screen changes when absorbing walls are present. 

7.1 A particle between two absorbing walls 

First, we consider a particle with symmetric absorbing walls at x = 0, a and zero 
potential. We assume that A_a = Aq. The wave function is given by 

^(a;,t)=|:A„exp|-^— ^tjsm-o;. 

It was shown in that for a particle with a single energy level the wave function decays 
at an exponential rate proportional to the energy. However, if there are more than just 
one level, the exponent contains beats. For example, for a two level system with real 
coefficients, we obtain the survival probability 

(7.1) 




The strongest beats occur for A; = 2, n = 1 with frequency ijJ\^2 = 3fnT'^/2ma'^. Setting 



Ai = A2 = Jl/2a and introducing the dimensionless time r = uji^2t^ we find that the 
survival probability is 

S{t) =exp|-^(5r-2sinr)|. 

The function 

log^(t) = -^(5r-2sinr) 



is qualitatively similar to that given in [15 



If absorption in windows of energy is assumed, the discounted wave function (|5.5|) is 



*(a;,t) = EA^expj-^— ^tjexpj ^^—^^Y^-^- (7-2) 
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A - — 



Writing 



such that X) kn| — I5 "W6 obtain that the survival probabihty (|5.6| ) is 



'' laJ^/inV 



Ssit) = Y^ |a„| exp 



-t . 



n=l 



ma'' 



(7.3) 



In the particular case that there is only one energy level in the box, both results (|7.1|) 
and ( ]7.cJ| ) should be the same. In this case, we obtain 



A = ira. 



7.2 A wave packet incident on an absorbing wall 



Next, we consider a Gaussian-like wave packet of free particles traveling toward an 
absorbing wall at a; = with positive mean velocity Ajq. The wave function is given by 



^(x,t) 



2a 



2\ 1/4 



n 






X 



(7.4) 



x + xo-^t) I I [x-xo + ^tj 
exp ^ — ^ ^ , 2iht + ^^0 {x + xo)} - exp { -^ , . ^,t^ ih {x - xq) 



a2 + ^^ 



a2 + ^^ 



where 



A 



1 — exp 



2xq 



exp 



2l2 



-a A;, 



and v? is a phase (see p, p.64]). It follows that 

2 



Xh 

rrni 



d_ 
dx 



^-(0,^) 



V2Xham'^ [Axl + a^fcg) 
(v/?)^ A Ua^m? + 4^2/^^ 



(7.5) 



exp —2a 



2 2 



2-^0 



+ U^klt^ — 2komxoth\ 



a^m? + Af^n' 



] 



hence 



Thus 



d_ 
dx 



^(0,t) 



dt < 00. 



R = Sioo) = lim [1 - Pit)] > 0, 



(7.6) 
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that is, the wave packet is only partially absorbed. This means that the "reflected" wave 
consists of trajectories that turned around before propagating past the absorbing wall. 
The absorption occurs when the trajectory propagates into the medium inside the wall. 
The discount of the wave function occurs when the packet is at the wall, as can be seen 
from eqs. ( |7.5| ) and ( |7.6| ). Thus R plays the role of a reflection coefficient. This is neither 
the usual reflection coefficient for a finite potential barrier nor that for an infinite barrier. 
The refiection coefficient is a function of the packet group velocity, ko, the width of the 
packet, a, its initial distance from the absorbing screen xq, the parameter A, and its mass 
m. This dependence is experimentally measurable. 

If absorption in energy windows is assumed, the discounted wave function, as given in 
eq. (|5.5|) , is calculated as follows. The modes are the functions 

ipk{x) = 
and the energies Ek are 




2 7.2 



Eh- 



h'k 



2m 

The discounted wave function is given by eq.( |5.5|) with summation replaced by integration 
with respect to k over the entire real line. The coefficient ak is the sine transform of the 
initial wave function ( [7.4| ) with t = 0, 



ak 




ikaP+ikoa^\ | ( f rf. —iko?+ikoa? \ 

J-O 2/1 I V 2 / 



e-P -^ exp -exp 



and the decay rate of the fc-th mode, as given in eqs. (|5.1| ), ( |5.5|) , is 

C I |2 i-2T 2 

tk \ak\ Ti k^ 



h 2inh 



ahe V27r r ky] { k^o?' 

exp < > exp 



cosh \kki^Q^\ — cos (2xoA;) 



4mA " [ 2 J " [ 2 ^ 
According to Section 5.1, the survival probability is the same as in the previous case. 

7.3 The slit experiment with lateral absorbing boundaries 

Here, we propose a simple device for performing a measurement of times of arrival of 
particles at an absorbing boundary in one dimension, as well as that of the unidirectional 
current at the absorbing boundaries. In addition, the proposed device demonstrates the 
effect of additional absorbing boundaries on the slit experiment. The proposed measure- 
ment can discriminate between the various modes of absorption described above. 
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Consider the setup of the sht experiment enclosed between two parallel absorbing 
walls, symmetric with respect to the slit and perpendicular to the planes of the screen 
and the slit. For example, the walls can be made of photographic plates. Particles are 
given a constant initial velocity, Vx-, in the x direction (perpendicular to the planes of the 
screen and slit), within the constraints of uncertainty. The time a particle leaves the slit 
is also measurable within the constraints of uncertainty. 

The initial packet is Gaussian in the x direction and is uniform inside the slit (in the 
y direction). This means that the initial velocities in the y direction have the density 



ii{k) 



sin ^k 



2^ 



(7.7) 



The plane of the slit is x = xq, the plane of the screen is x = 0, the slit is the interval 
— 7r/2 < y < 7r/2. The absorbing planes are y = ±yo with yo > Tr/2. In this setup the 
motion of the particles in the x direction is independent of that in the y direction. The 
latter is the object of the proposed experiment. 

Particles that hit the planes y = ±|/o leave traces at points xi,X2, ■■.,xn. These 
distances are proportional (within the constraints of uncertainty) to the times of arrival 
at the absorbing walls of particles that start out in the interval —tt/2 < y < 7r/2 with 
initial velocities distributed as in eg. ([7. 71). The histogram obtained from these points, on 
an axis normalized with the velocity Vx, is that of the times of arrival of one dimensional 
particles moving on the y axis. 

The wave function for this configuration is given by 

^ix,y,t) = ij^{x,t)^lj^{y,t), 



,2 . 



where \dtlj^{0,t)/dx\ is the same as given in eq.( |7.5| ) and 



,2/ .N 'F^ 2 nir^ . rnr 

V [y, t) = 1^ —JU^ cos -J— sm —y exp <; -- 



i ra^TT^ 



n=l 



nvr 



3/2 



4^/0 



1/0 



^ Vq 



According to our theory, the pattern on each lateral wall is given by 



J{x, ±yo 



J{x,±yo,t)dt 



dip^{x,t) 



dx 



V^'(±yo,t) 



dt. 



(7. 



Next, we consider the pattern observed on the wall x = 0, — 7r/2 < y < 7i/2. The 



histogram of the traces of the particles on the screen at x = is given by 

2 



Jiy) 



dtlj\0,t)/dx\ 



E 

n=l 



riTT 



3/2 



niT nir 

cos sm — y exp < — 

4yo Vo I 



i n^TT^ 



yi 



dt. (7.9) 
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If the velocities in the x direction are concentrated around Vx-, the histogram will be 
approximately 



J{y) 



J{y,l)= aV'^(0,t)/ax ^lj\y,f) 



d'ilj\Q,f)/dx\ 



E 

n=l 



2 nvr nvr 

— rrr COS - — sm — y exp 

rni^y^ %o yo 






This is not the same as the expression obtained in eg .( 16.81) . The difference is due to 
the effect of the lateral absorbing boundaries. Thus, according to our theory, absorbing 
boundaries cannot be ignored, as usually done in quantum mechanics. 

If absorption in energy windows is assumed, the function %lj^{x,t) is the same as cal- 
culated in Section 7.1 above, 



W [y, t> = }^ zrz^ cos -J— sm —y exp 



i n^vr^ 



An vr 



2^2 



n=l 



niT 



3/2 



4z/o 



Vo 



^ Vo 



yl 



The patterns on the absorbing walls are given by eqs. (|7.8|) and (|7.9|) . 
8. Discussion and summary 



We observe that according to eq. ( [4.5| ) there is conservation of probability: the prob- 
ability of the absorbed Feynman trajectories and that of the surviving trajectories sum 
to 1. This is the result of our postulate that the Feynman integrals over the two classes 
of trajectories have disjoint supports, that is, the absorbed trajectories never return to 
the interval [a, 6]. This conservation of probability persists for all times. This result is 
different from that obtained in decoherent state theory ||^ . 

The survival probability S{t) is the probability obtained by repeating the experiment 
of observing the particle between the absorbing boundaries at time t and constructing 
a histogram of the number of times the particle is observed. This probability is not 
a quantum mechanical quantity in the sense that it is not the integral of the squared 
modulus of a probability amplitude defined by Schrodinger's equation. 

The wave function of the particle at time t, given that it has not been absorbed by 
that time, is the wave function ^(x, t) defined by Schrodinger's equation inside the given 
domain with reflecting boundary conditions. That is, if the absorbing boundary represents 
a detector, the wave function of the particle is ip{x, t) as long as the particle has not been 
detected. 

The concept of absorption of particles, energy, momentum, and so on, is an aspect 
of time irreversible processes in quantum mechanics. Examples of irreversible processes 
in quantum mechanics are the notion of collapse of the wave function that is caused by 
measurement (this will discuses at |]TB| ), the decay of a state of a particle, a particle that 
enters into a bulk and loses its momentum due to interactions with the particles of the 
bulk, and so on. This paper is an attempt to construct a formalism for the description of 
such phenomena. 
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The process of absorption of quantum particles can be illustrated by a packet fo iden- 
tical non-interacting particles that hit an absorbing surface, for example, a photographic 
plate. At the moment a particle hits the plate it is absorbed in the sense that its wave 
function no longer evolves according to the Hamiltonian of the particles that have not 
been absorbed so far. Thus the absorbing surface separates the particles of the packet 
into two sets, those that have reached the surface by a given time and those that have 
not. The interference between the Feynman integral over the trajectories of the absorbed 
particles and that over the trajectories of the surviving ones vanishes. The probability of 
the absorbed particles is discounted from tha total probability of the packet. This is also 
the case for any other detector that absorbs particles. The surviving trajectories, those 
that have not reached the absorbing boundaries so far, give rise to a reflected wave, as 
if the absorbing boundary were an infinite potential wall. This fact becomes apparent 
not from a solution of a wave equation, but rather from the calculation of the Feynman 
integral over a class of restricted trajectories. 

Our derivation does not start with a Hamiltonian, but rather with an action of tra- 
jectories in a restricted class. 

Quantum mechanics without absorption is recovered from our formalism when the 
absorbing boundaries are moved to infinity or when the absorption constant A vanishes. 

The examples demonstrate the expected phenomenon that particles that reach the ab- 
sorbing boundary are partially reflected and partially absorbed. In either case the decay 
pattern of the wave function seems to be new. 

Acknowledgment: The authors wish to express their gratitude to Y. Aharonov and L. 
Horowitz for useful discussions. 
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